Coherent perfect absorbers for transient, periodic or chaotic optical fields: 
time-reversed lasers beyond threshold 
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Recent works [Y.D. Chong ei ai, Phys. Rev. Lett. 105, 053901 (2010); W. Wan et al, Science 
331, 889 (2011)] have shown that the time-reversed process of lasing at threshold realizes a coherent 
perfect absorber (CPA). In a CPA, a lossy medium in an optical cavity with a specific degree of 
dissipation, equal in modulus to the gain of the lasing medium, can perfectly absorb coherent optical 
waves at discrete frequencies that are the time-reversed counterpart of the lasing modes. Here the 
concepts of time-reversal of lasing and CPA are extended for optical radiation emitted by a laser 
operated in an arbitrary regime, i.e. for transient, chaotic or periodic coherent optical fields. We 
prove that any electromagnetic signal E{t) generated by a laser system S operated in an arbitrary 
regime can be perfectly absorbed by a CPA device S' which is simply realized by placing inside S 
a broadband linear absorber (attenuator) of appropriate transmittance. As examples, we discuss 
CPA devices that perfectly absorb a chaotic laser signal and a frequency-modulated optical wave. 

PACS numbers: 42.55.Ah, 42.65.Sf 



I. INTRODUCTION 



A laser oscillator is a device that self-organizes to emit 
a narrow-band coherent electromagnetic radiation when 
the pumping level exceeds a threshold value. Above the 
first lasing threshold, lasers are generally complex non- 
linear systems [Ij and may emit different kinds of output 
signals E{t), ranging from the ideal monochromatic wave 
under single-mode continuous-wave operation, to irregu- 
lar or chaotic signals when multi-mode or other kinds 
of instabilities set in [1, 2j, or giant pulses or periodic 
trains of ultrashort optical pulses when operated in the 
Q-switching or mode-locking regimes [3l |4j. Generally 
speaking, a "time-reversed" laser refers to a device that, 
rather than emitting the signal E{t) propagating outgo- 
ing from the laser cavity, it is capable of perfectly ab- 
sorbing the same signal E{t) that propagates backward 
into the cavity, without any reflection [5j. Such a de- 
vice thus realizes a coherent perfect absorber (CPA) for 
the signal E{t). For a continuos-wave field, it is known 
that a dissipative resonator can perfectly absorb inci- 
dent light at the condition of critical coupling [B]. Re- 
cently, Douglas Stone and coworkers have suggested and 
proven rather generally the possibility to realize a CPA 
exploiting the time-reversed process of lasing at thresh- 
old. In this case, the laser system behaves as a linear 
one, and the time-reversal process is realized simply by 
the replacement e(r) — ^ e*(r) for the complex dielectric 
constant e of the medium. In this way, any linear ab- 
sorbing medium of arbitrary shape behaves as a CPA at 
some discrete frequencies under appropriate coherent il- 
lumination (the time reverse of the output lasing modes) 
and provided that a precise amount of dissipation in the 
medium occurs (equal in modulus to the threshold gains 
for lasing). The CPA idea has received a great inter- 
est and stimulated theoretical and experimental stud- 
ies along different lines [71116). In particular, in Ref.[S] 



an experimental demonstration of interferometric control 
of the absorption based on CPA was reported using a 
thin slice of silicon illuminated by two beams, whereas 
in Refs.[71[n] a laser-absorber device, which can operate 
as a CPA and as a laser simultaneously, has been sug- 
gested combining the CPA and T'T-symmetry concepts. 
Time-reversal of other optical instabilities, such as time- 
reversal of optical parametric oscillation, have been pro- 
posed to realize a multicolor CPA in Ref.|llj. Also, ex- 
tension of the concept of time-reversed laser and CPA to 
the spaser and plasmonic nanostructures have been sug- 
gested in Refs. [131 E] • Here the energy of the incoming 
wave is fully transferred into surface plasmon oscillations 
and evanescent electromagnetic fields. CPA-based de- 
vices may have potential applications to the realization 
of a new class of absorptive interferometers and nanosen- 
sors. In all previous studies, time-reversal of lasing has 
been limited to consider either a laser at threshold or 
above threshold in steady-state operation. However, as 
previously mentioned, a laser can operate in rather com- 
plex or transient regimes, which are highly nonlinear. A 
major and foundational open question is whether there 
exists a device that realizes the time-reversal of a laser 
operating in any (generally highly-nonlinear) regime, i.e. 
capable of perfectly absorbs the field E{t) emitted by a 
laser operating in any regime. If yes, how can we realize 
(at least in principle) such an "anti-laser" device? 
It is the aim of this work to answer to such two ma- 
jor questions. By considering an optical system with a 
single input /output channel and in the plane- wave ap- 
proximation, we will prove rather generally that for any 
electromagnetic signal E{t) generated by a laser system 
S operated in an arbitrary regime, i.e. emitting a tran- 
sient, irregular, chaotic, or periodic signal, one can al- 
ways construct a CPA device S' that perfectly absorbs 
the field E{t) emitted by S, and that a possible simple 
realization of this system is obtained by placing inside S, 
nearby the output coupler, a broadband linear absorber 



(attenuator) of appropriate transmittance. As examples, 
we discuss two CPA devices that perfectly absorb the 
former a chaotic signal emitted by a single-mode laser 
operated in the Lorenz-Haken instability regime, the lat- 
ter a frequency-modulated (FM) optical wave emitted by 
a multimodc FM-operated laser. 



II. BASIC CPA IDEA 

Let us consider a rather general ring-cavity laser sys- 
tem S of length £ with a single output coupler, consist- 
ing of a lossless beam splitter BS, as shown in Fig. 1(a). 
The laser cavity contains a gain medium and possibly 
other optical elements or devices, such as saturable ab- 
sorbers, amplitude or frequency modulators, etc. de- 
pending upon the operating regime of the laser. The 
electric field £{z,t) circulating inside the cavity can be 
written as £{z, t) — A(z, t) exp(ikz—ujt), where cj is a ref- 
erence frequency, k — uj/c is the wave number in vacuum, 
z is the longitudinal spatial coordinate along the ring, 
and A{z,t) is a slowly-varying envelope. Without loss of 
generality, a; is chosen to be a resonance frequency of the 
empty cavity, so that kC is an integer multiple of 2tt. Af- 
ter one cavity transit, the envelope (/)(i) = A{z = L^ ,t) 
at the plane z — C~ can be formally written as 



4>{t)^P{^{t~Tj,),t), 



(1) 



where -0(0 — ^(^ = 0+,i) is the field envelope at the 
z = plane, T/j = C/c is the cavity transit time, and 
P is a functional of il}{t — Tr) that describes the field 
evolution over one cavity round trip. The specific form 
of P(V': t) depends on the optical elements and devices 
put inside the cavity and in general can not be given 
explicitly, requiring to solve coupled equations that ac- 
count for e.g. polarization and population dynamics in 
the active medium or in saturable absorbers possibly 
present inside the laser cavity with the appropriate ini- 
tial conditions for the population and polarization vari- 
ables (see, for instance, [T7]). Exphcit forms of the op- 
erator P can be given is some cases, for example, for 
mode-locked lasers [18]. However, for our purposes we 
do not need here to explicitly specify the form of P. 
The time delay Tr in Eq.(l) is introduced in such a way 
that for an empty cavity P is the identity operator, i.e. 
P{ip{t-TR),t) ^ip{t~TR). The evolution of ■0(i) at suc- 
cessive cavity transits and the corresponding output laser 
field E{t) are simply obtained by imposing the scatter- 
ing relation between input/output channels at the beam 
splitter BS. For a lossless BS, the scattering matrix is 
unitary, and one can write (see the inset of Fig.l) 
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R 
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(2) 



where T and R = 1 — T are the BS transmittance and 
reflectance, respectively, that are assumed to be spec- 
trally flat. For the laser system S without an injected 



signal, one has f{t) — 0, 4){t) = {1 / \/ R)'ip{t) and thus 
one obtains the following equation for ip{t) 



4>{t)^VRPmt-TR),t). 



(3) 



The signal E{t) emitted by the laser system in then given 

by 



E{t) 




(4) 



Once the initial field distribution ^p{t) is assigned in the 
interval < t < Tr, Eq.(3) can be used to determine the 
evolution of ipit) at successive transits in the ring. To 
clarify this point, let us introduce a local time variable 
T, with < T < Tfj, and let us set t = nTu + r, where 
n = 0,1,2,3,... is the round trip number in the ring. 
After setting ■(/''"'(t) = ^{t = tlTr + r), Eq.(3) can be 
written as 



^(")(r) = yPP (V^""'^ (t), r + uTr^ 



(5) 



For an assigned initial condition '(/'('^'(t), the map 
(5) enables to determine recursively the field distri- 
butions iP^^\t) at successive round trips. The func- 
tional P (?/''^"^^)(r),T -I- tiTr) is generally computed in 
the frame of an initial-boundary value (Goursat) prob- 
lem of coupled differential equations (see, for instance, 
[17j : see also the example discussed in the Appendix). 
The temporal evolution of the emitted field ^p{t) depends 
on the operating regime of the laser. For example, it 
could describe a periodic train of ultrashort pulses when 
the laser is operated in the mode-locking regime; in this 
case the self-consistent equation (3) has the form of a 
differential-delayed equation, the so-called master equa- 
tion of mode- locking |21 HB] . But it could also describe 
a rather irregular waveform, such as for a laser operated 
in the chaotic regime (see, for instance, O US]). The 
main question is now: is there an optical system S' that 
can perfectly absorb the optical field E{t) emitted by the 
laser system S, i.e. can we in some way time reverse a 
laser that is operated in an arbitrary regime? 
The answer is affirmative, and a possible realization of 
the system S' is depicted in Fig. 1(b). S' is basically 
obtained from S by just inserting, nearby the output 
beam splitter BS, a broadband linear absorber (atten- 
uator) with a transmittance Ta given by 



Ta = R'. 



(6) 



In this way, it can be shown rather generally that S' 
behaves as a CPA when the signal fi{t) injected into the 
cavity is given by 



(7) 



h{t) = ^RE{t) = VTi,{t) 



where E{t) is given by Eq.(4) and ijj{t) is the solution 
to Eq.(3) with the appropriate initial conditions. Such 
an injected signal is basically the output of the laser S, 
attenuated by the factor \/P [2D]. 
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FIG. 1. (Color online) Schematic of (a) a laser system S, 
made of a ring-cavity of length £. with a single output cou- 
pler (lossless beam splitter BS, with transmittance T and re- 
flectance R — 1 — r), and (b) corresponding CPA system S'. 
Note that S' is simply obtained from S by the insertion of a 
linear absorber (attenuator) nearby the output coupler with 
transmittance Ta = Fl? . The inset in (a) shows the scattering 
relations of the field amplitudes at the lossless BS. 



To prove that S' time reverses the laser system S when 
Eqs.(6) and (7) are satisfied, let us inject into S' via the 
beam splitter port the signal fi{t). Let us then indicate 
by f 1 (z, t) = Ai{z, t) exp{ikz — itut) the electric field that 
is established inside the cavity of the system S' after 
signal injection, and by Ei(t) the corresponding output 
field [see Fig. 1(b)]. In this case one can obviously write 
[compare with Eq.(l)] 



Mt)^VTaP{Mt-TR),t) 



(8) 



where we have set ipi{t) = Ai{z = 0+,t) and (/'i(t) = 
Ai(z = C~,t). In fact, the cavity round-trip operator 
A of the system S' is simply obtained by cascading the 
propagator P of S with the transmission amplitude ^/Ta 
of the added linear absorber, i.e. Pi = y^TaP. On the 
other hand, the scattering relation of the field amplitudes 
at the beam splitter BS in the system S' relates V'i(Oj 
(j)i{t) and the injected signal /i(i) according to 



(9) 



Mt) = VTfi{t) + VRMt) 



From Eq.(9) one obtains (j)i{t) = {l/VR)^pi{t) - 
\/TjRfi{t), which after substitution into Eq.(8) yields 
the following equation for the intracavity field "01 {t) that 
is established inside the cavity at the plane z = 0+ 

Mt) = Vfh{t) + y/RTaP{Mt~TR),t) (10) 



The field Ei{t) leaving S' is then obtained from the scat- 
tering relation at the beam splitter and reads explicitly 
Ei{t) = -y/Rfi{t) + VTcj)i{t), which by means of Eq.(9) 
can be written as 



Ei{t) = 



IT . ,.. hit) 



(11) 



Similarly to the case of the laser system S in the ab- 
sence of the injected signal, Eq.(lO) can be written as a 
recursive relation 

V;|")(r) = ^^/["'M + ^/RTaPi4''-'\T),T + uTr) 

(12) 
where we have set f/'i (t) — i^iir + uTr) and /i (t) — 
/i {nTR + t) . Once the initial field distribution ^J (r) is 
assigned, the map (12) enables one to determine the evo- 
lution of the intracavity field -01 (0 ^-t successive round- 
trips. Let us now assume that the two systems S and 
S' are initially prepared in the same state, i.e. that 
■0j^ '(t) = ■(/'("^(r) and that atomic variables (population 
inversion and polarization) have the same initial values 
in the active medium (for example their equilibrium val- 
ues). In this case, one can readily show using Eq.(5) that, 
if the injected signal /i(i) is chosen according to 



hit) 



rr, 



-V'(t) 



(13) 



the recursive relation (12) admits of the solution 

Vi {'^) — V'^"'' (''") for any n. This implies that the intra- 
cavity fields V'i(i) and -0(0 established in the two systems 
S and S' are the same. Note that this result holds even 
if the map (12) admits of different attractors (such as 
for a bistable system), provided that the initial condi- 
tions in S' belong to the basin of attraction of the CPA 
solution ipi [t) = ?/;("^(t). This ensures that asymp- 
totically one has 05;" (r) -^ ■0^")(r) as n ^ cx), even 
though 0]^" (0) 7^ '0''"HO)- But for some other initial 
conditions it might happen that "0}" (r) is not attracted 
toward 0'^"^(r), and the CPA scheme fails. However, like 
in any nonlinear dynamical system showing different at- 
tractors, one can switch the path from one stable attrac- 
tor to another one by introducing large perturbations of 
the system parameters (e.g. by transiently changing the 
cavity losses or the gain in the medium) . As discussed in 
Sec. III. A with reference to a specific example , the result 
holds for transient or chaotic regimes as well, which are 
very sensitive to the initial state of the system. 
The field Ei (t) emitted by the system S' is obtained after 
substitution of Eq.(13) into Eq.(ll). For ■0i(i) = ip{t), 
one obtains 



Ei{t) 



/Ta-R 



^{t) 



(14) 



which vanishes if the transmittance T^ is chosen to sat- 
isfy Eq.(6). Correspondingly, the injected field /i(i), as 



obtained from Eq.(13) after setting Ta = R^, is given by 
Eq.(7). Hence, provided that Eqs.(6) and (7) are met 
and the two systems S and S' are initially prepared in 
the same state, S' behaves as a CPA device for the signal 
emitted by S, regardless of its operational regime. 



III. EXAMPLES OF CPA DEVICES 

In this section we discuss two examples of CPA de- 
vices corresponding to a laser oscillator S operating in 
two nontrivial regimes. The former example is a CPA de- 
vice for a chaotic optical signal emitted by a single-mode 
homogeneously-broadened laser operated in the chaotic 
(Lorenz-Haken instability) regime [UHIIIF ; the latter ex- 
ample is a CPA device for a frequency-modulated (FM) 
optical signal emitted by a FM-operated laser with an 
intracavity phase modulator [3l [22j [23] . 



A. CPA for a chaotic optical field 

The first example we would like to discuss is a CPA 
device that perfectly absorbs a chaotic optical field emit- 
ted by a single-mode homogeneously-broadened laser op- 
erated in the chaotic (Lorenz-Haken instability) regime 
[HEIIIS]- In this case, the laser system S just contains 
a homogeneously-broadened two-level medium of length 
I with a small-signal gain coefficient per unit length g, 
which is provided by population inversion in the medium. 
In the single longitudinal mode and uniform field approx- 
imations, the intracavity field A{z,t) is assumed to be 
almost uniform along the ring (i.e. almost constant for 
< z < £) and slowly-varying in time over one cav- 
ity round-trip time j5]. In this case, the map (5) can 
be effectively replaced by a set of three coupled differen- 
tial equations describing the evolution of the intracavity 
field, polarization and population inversion in the two- 
level medium, the so-called Lorenz-Haken model of the 
single-mode homogeneously-broadened laser (see the Ap- 
pendix for technical details). These equations read |2j 



dtp 

lit 
dA 

lit 
dn 

~di 



-n(i> + 2CA) 

-7j_[(l + iA)A-|-V^n 
1 



(15) 



-711 



n-l 



(V^A*+V^*A) 



where: '(p{t), A(i) and n{t) are the intracavity electric 
field, polarization and population inversion, respectively, 
normalized as in Refs.L2j J^2J; 711 and 7_l are the popula- 
tion and dipole decay rates, respectively; 



Tc 
2£ 



(16) 



is the cavity decay rate; C = gl/T; and A = (wq — 
^)/"f± is the normalized detuning parameter between the 
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FIG. 2. (Color online) (a) Evolution of the intracavity field 
tp{t) in the laser system S as obtained by numerical solution 
of the Lorenz-Haken equations (15) for the initial condition 
(V'(O) = 0.001, A(0) = 0,n(0) = 1) (solid curve), and (b) 
corresponding projection of the phase-space trajectory in the 
('i/',n) plane. Parameter values are n/'yi_ ~ 4, 7||/7± = 0.5, 
2C = 40 and A = 0. In (a) the dotted curve shows the 
evolution of ip{t) as obtained by the sHghtly different initial 
condition {ip{0) = 0.0011, A(0) = 0,n(0) = 1). Note the 
strong sensitivity on the initial condition, which is a rather 
general feature of a chaotic attractor. In (c) we show the nu- 
merically computed evolution of the output field Ei (t) , nor- 
malized to VT (i.e. the behavior of ipii^') ~ ''Pi^))j in the 
CPA system S' with an injected field fi{t) and attenuator 
transmittance Ta satisfying Eqs.(6) and (7), for the initial 
conditions {ipi{0) = 0.0001, Ai(0) = 0,ni(0) = 1) (curve 1), 
and {ipi{0) = 0.005, Ai(0) = 0,ni(0) = 1) (curve 2). The 
dotted curve in (c) shows, for comparison, the behavior of 
the injected field fi{t), normalized to VT [i.e. the behavior 
of 4'i't)]- (d) Same as (c), but for a large deviation of initial 
conditions {Tpi{0) = 5, Ai(0) = 0,ni(0) = 1). 



cavity resonance frequency a; and the atomic transition 
frequency cjq. In particular, at resonance A = (i.e. for 
uj — ujq), a and ip can be taken to be real-valued, and 
Eqs.(15) are analogous to the Lorenz model, developed 
for convective instabilities in hydrodynamics [UIJ]. In the 
following, we will mainly limit to consider the limit A — 
0. In this case, laser threshold is attained for 2C > 1, at 



which the non-lasing trivial solution {tjj = 0,A = 0,n = 
1) to Eqs.(15) becomes unstable. Above threshold, the 
laser equations admit of a steady-state solution, with two 
possible phases for ^/j and A (0 and tt), namely 



1 



V'i = ±V2C - 1, A± = T^V2^- 



1, n± = 



1 
2C' 



(17) 

Such solutions may undergo a subcritical Hopf instabil- 
ity in the bad cavity limit (k > 711 -I- 7j_) and for large 
small-signal gains [D 12] ■ In such a regime chaotic self- 
pulsations can be observed. Here, the point (ipit), A(i), 
n{t)) in phase space never settles down but continually 
makes excursions about one of the two laser solutions 
(17) with what appear to be random jumps from cir- 
cling one fixed point to circling the other. The strange 
set is known as the Lorenz attractor, and the output 
laser field E{t) ~ \/Tij}{t) turns out to be strongly sen- 
sitive to the initial conditions [H H]. As an example, 
in Fig. 2(a) we show the evolution of the intracavity laser 
field ^{t) as obtained by numerical solution of the Lorenz- 
Haken equations (15) for parameter values k/7_l — 4, 
7||/7-L = 0.5, 2C = 40, A = and for the initial con- 
dition ?/.(0) = 0.001, A(0) = 0, n{Q) = 1. The cor- 
responding projection of the phase-space trajectory in 
the {ip,n) plane is depicted in Fig. 2(b), which shows the 
characteristic Lorenz attractor of the laser dynamics in 
the chaotic regime. It should be noted that in such a 
regime the established intracavity laser field (and hence 
the laser output field) is strongly sensitive to the initial 
conditions. For example, the dotted curve in Fig. 2(a) 
shows the intracavity field that one would observe for 
the slightly changed initial conditions ip{0) = 0.0011, 
A(0) = 0, n(0) = 1. 

Let us now consider the CPA device S' associated 
to S, as discussed in Sec.IL In this case, in the single- 
longitudinal mode and mean-field approximations, the 
dynamical evolution of the intracavity field V'l (t) is gov- 
erned by a set of coupled equations sinfilar to Eqs.(15), in 
which the first equation is modified to take into account 
the effects of the injected field fi{t) and of the intracav- 
ity absorber with transmittance Ta . One obtains (see the 
Appendix for technical details) 



dt 
dAi 
"dt 

dn 



= -k(VJi + 2CAl) - Kltpl + K 

= -7_L[(l + iA)Ai+V'ini) 
1 



2/1 (i) 



(18) 
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n-1- -(M* +^*A) 



where 



Ki 



c(l 



/t: 



,(1 



/t:, 



c 



T 



(19) 



The output field Ei (t) from the CPA device is then given 
by Eq.(ll). According to the general analysis presented 
in Sec. II, if the injected field is chosen to satisfy the 



condition (13), it can be readily shown that the solu- 
tion to Eqs.(18) is given by ipiit) = ip{t), Ai(t) = A(t), 
ni{t) — n{t) provided that the same initial conditions 
V'i(O) = VIO), Ai(0) = A(0), ni(0) = n(0) are assumed 
for the field and atomic variables in the two systems S 
and S'. Furthermore, if the transmittance Ta of the ab- 
sorber in S' is tuned to satisfy the condition (6), the 
output field Ei{t) emitted by S' vanishes, i.e. S' behaves 
as a perfect CPA for the chaotic optical field ij{t). Note 
that for the perfect CPA device one has ki = 2k and the 
dynamical system (18) reads 



chpi 
dt 
dAi 
"dt 
dni 
~dt 



-k(i/;i -I- 2CAi) - 2k[V'i - V(0] 
-7_L[(l + iA)Ai+V'ini) 



(20) 
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ni-l--(V'iA*+ViAi) 



where '4){t) is the solution to Eqs.(15) with the assigned 
initial conditions. 

It should be noted that in a practical case the laser 
system S and the CPA system S' can not be exactly 
prepared in the same state. In this case, the solu- 
tion (■0i,Ai,ni) to Eqs.(20) deviates from the solution 
(i/). A, n) of Eqs.(15), just because the initial conditions of 
field and/or atomic variables are different. Correspond- 
ingly, the output field Ei{t) emitted by the CPA device, 
given by 



Ei{t) = \l j^[i^i{t) 



^(t)]^VT[^i(i)-^(i)], (21) 



would not vanish, i.e. perfect absorption would be lost. 
Such a circumstance rises the question whether the CPA 
idea presented in Sec. II is actually of physical relevance. 
This objection is especially a serious one whenever the 
system dynamics is strongly sensitive to the initial condi- 
tions, such as for a chaotic laser, or if the system S' with 
injected signal may show two (or more) stable attrac- 
tors, i.e. in the presence of bistability or multistability 
|21) . In this case, an initial condition of S' sufficiently far 
from that of S can bring the dynamical system (20) into 
a different attractor than {tp{t),A{t),n{t)) [H]. Strong 
deviations of the initial conditions may occur, for exam- 
ple, whenever the pump parameter 2C is larger than 3 
and system S' thus self-oscillates, i.e. it is above laser 
threshold (in the absence of the injected signal) in spite 
of the attenuator put in the cavity. However, provided 
that the basin of attraction of {ip{t),A{t),n{t)) for S' 
is sufficiently wide, a large perturbation can switch the 
path of the dynamical system toward the "right" CPA 
attractor {ip{t), A{t), n{t)). Of course the basin of attrac- 
tion of {ip{t),A{t),n{t)) for S', as well as the kind and 
strength of the large perturbations requested to switch 
the nonlinear dynamics of Eqs.(20) from one attractor 
to another one, should be considered on a case-by-case 
basis. Let us focus here our attention to the Lorenz- 
Haken model of laser chaos. In this regime it is known 



that even a small change in the value of '0(0) niay deeply 
modify the output waveform, as shown in Fig. 2(a). How- 
ever, numerical results show that the basin of attrac- 
tion of (0(t), A(i),n(t)) for Eqs.(20) is quite broad, i.e. 
asymptotically one has ipiit) ~^ ''Pit) for ^ quite broad 
range of initial conditions (01(0), Ai(0), ni(0)) around 
(0(0), A(0), n{0)). The reason thereof is that the dynami- 
cal system (20) [contrary to the Lorenz-Haken dynamical 
system (15)] is driven by the external field 0(i), which 
forces ■0i(*) toward 0(i) after an initial transient. This is 
clearly shown in Fig. 2(c), which depicts the output field 
Ei{t), normalized to a/T [i.e. the difference 0i(i) — i^{t); 
see Eq.(21)] emitted by the CPA device S', as obtained 
by numerical integration of Eqs.(20), for the two ini- 
tial conditions (V'i(O) = 0.0001, Ai(0) = 0, ni(0) = 1) 
[curve 1 in Fig.2(c)] and (0i(O) = 0.005, Ai(0) = 0, 
ni(0) = 1) [curve 2 in Fig. 2(c)], which appreciably differ 
from 0(0) = 0.001, A(0) = 0, 7i(0) = 1. For comparison, 
in the figure we show (but on a different scale) the field 
fi{t) — vTil>{t) injected into the cavity, normalized to 
^/T (i.e. the behavior of Tp{t)), where ijj{t) is the field 
of S for the initial condition (0(0) = 0.001, A(0) = 0, 
n(0) = 1) [i.e. the solid curve in Fig. 2(a)]. Note that, 
since the initial conditions in S' and S are different, the 
CPA output field does not exactly vanish in an initial 
transient [see curves 1 and 2 in Fig. 2(c)], however it re- 
mains much smaller than the injected field even though 
the initial conditions in the two systems are appreciably 
different. Most importantly, after an initial transient the 
output field vanishes, and perfect absorption is attained. 
Such a behavior is observed even for large deviations of 
initial conditions, as shown in Fig. 2(d) as an example. 



B. CPA for a frequency-modulated optical field 

As a second example, we discuss the realization of 
a CPA device for an optical frequency-modulated sig- 
nal E[t), i.e. an optical field with a constant intensity 
but with a sinusoidally-modulated optical phase. A co- 
herent FM optical signal E{t) is approximately gener- 
ated by placing a phase modulator inside a laser device 
and sinusoidally-driving it asynchronously as compared 
to the cavity round-trip time. In this way the laser oper- 
ates in the so-called frequency-modulation (FM) regime 
[31 [221 [23] • The FM regime is a multimode regime 
in which several cavity axial modes are excited out of 
resonance with steady-state amplitudes but with time- 
varying phases, resulting in a constant intensity laser field 
but with a carrier laser frequency which is sinusoidally 
swept at the frequency impressed by the phase modulator 
[3] . As shown in Ref . [23^ , a pure FM signal at optical car- 
rier can be generated provided that the phase modulation 
is sufficiently asynchronous such that the bandwidth of 
the FM signal is much smaller that the gain bandwidth of 
the active medium. A schematic of a FM-operated laser 
is shown in Fig. 3(a). The optical cavity contains a gain 
medium and a phase modulator, which is sinusoidally- 



driven at a frequency ujm which is enough detuned from 
the cavity axial mode separation ujax = 27r/TR, where 
Tr is the cavity photon transit time [3] . For a slow gain 
medium and neglecting dispersive and finite gain band- 
width effects, the cavity round-trip operator reads [23 

P(i,0) = exp(g)exp[iAcos(w™i)]0(i- Tfl), (22) 

where g is the single-pass saturated gain in the active 
medium and A is the modulation depth impressed by 
the phase modulator. The self-consistent equation for 
the field Tp{t) [Eq.(3)] then reads 

0(i) = exp[.g - 7 + iAcos(w„t)]0(i - T^), (23) 



where 7 = — Invi? is the logarithmic loss of the cav- 
ity due to the output coupling. Equation (23) should 
be supplemented with a rate equation for the saturated 
gain g (see, for instance, [24j). After an initial transient, 
a steady-state operation is achieved, in which the sat- 
urated gain settles down to the stationary value .9 = 7 
[24|. Correspondingly, the intracavity field '0(i) is given by 
one of the "modes" of the phase- modulated cavity (gen- 
erally the one with the carrier fequency closest to the 
center of the gainline) . Such modes are simply obtained 
as eigenf unctions of Eq.(23) and reads explicitly [23] , [24 ] 

ip{t) — exp[ircos(cL'm^ + ¥') + ill^axt]- 



(24) 



(26) 



In Eq.(24), I = 0,±1,±2,... is the mode index, F is the 
effective modulation index, and ip a phase offset. The 
values of F and </? are obtained after substitution of the 
Ansatz (24) into Eq.(23) and equating the imaginary 
terms in the exponentials on the left and right hand sides 
of the equation so obtained. This yields 

r = TT-^, 1 — ^ (25) 

2sm.['n:ujm/'^ax) 

A 
cos if = — . 

^ 2r 

The output field, ~ 0(t)e*"*, is thus a pure FM signal, 
i.e. a frequency comb, with spectral lines at frequencies 
nujrn [n = 0,±1,±2, ...) around the carrier frequency 
Lj + lujax of amplitudes ^ J„(r), where J„ is the Bessel 
function of first kind and of order n, and w is a reference 
cavity resonance frequency. According to Eq.(25), the ef- 
fective modulation index F increases and diverges as the 
synchronous modulation condition Wm = oJax is attained. 
Indeed, in practical FM-operated lasers large modulation 
indices are usually obtained by tuning the modulation 
frequency LOm close to the cavity axial mode spacing ujax 
(see, for instance, [22 [2H])- However, the solution (24) 
provides an accurate approximation to the output field 
of the FM-operated laser provided that the modulation 
frequency Um remains sufficiently detuned far apart from 
LUax in such a way that the spectral extent ~ 2ra;m of the 
FM signal ^(i) is much smaller than the gain linewidth 
ujg, i.e. provided that the following condition 
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FIG. 3. (Color online) Schematic of (a) a FM-operated laser 
S, made of a ring-cavity containing a gain medium (gain co- 
efficient g) and a phase modulator, driven at a frequency ujm 
detuned from the frequency separation uuax of cavity axial 
modes, and (b) the corresponding CPA system S', which is 
obtained from S by replacing the gain medium with an ab- 
sorber with absorption coefficient —g. 



is satisfied. Indeed, as the synchronous modulation con- 
dition LOm = Wax IS approachcd, laser operation switches 
into a pulsed regime (FM mode-locking), with the gen- 
eration of a periodic train of short optical pulses [3l |26] . 

Let us now discuss the possibility to perfectly absorb 
the FM signal emitted by S. According to the analysis of 
Sec. II, a CPA system S' that perfectly absorbs the FM 
signal can be obtained from S by placing inside the cavity 
a linear absorber with a transmittance Ta — R^. Note 
that, as 5 = 7 = —\n\/R and since we neglected gain 
bandwidth and dispersion effects, the combined effect of 
the gain and absorber media in S' is equivalent to that 
of a single linear absorber placed in the cavity with an 
absorption coefficient —g, as shown in Fig. 3(b). In other 
words, for this special example we retrieve the simple rule 
of Ref. [F that the CPA system is obtained from the lasing 
system by replacing the gain mediunr with an absorber 
with an amount of dissipation exactly opposite to the 
amplification factor in the lasing medium at threshold. 
We stress, however, that the CPA system S' that realizes 
the time reverse of the lasing system S, obtained by the 
general procedure outlined in Sec. II, does not generally 
correspond to the replacement of the gain medium with 
an absorber, as discussed for instance in the example of 
Sec.III.A. 



IV. CONCLUSION 

In conclusion, in this work we have extended the con- 
cepts of time-reversal of lasing and CPA, recently intro- 
duced in Ref. [S] , to the case of optical radiation emitted 
by a laser operated in an arbitrary regime, i.e. for tran- 
sient, chaotic or periodic coherent optical fields. We have 
proven rather generally that any electromagnetic signal 
E{t) generated by a laser system S operated in an ar- 
bitrary (and generally highly-nonlinear) regime can be 
perfectly absorbed by a CPA device S' which is simply 
realized by placing inside S a broadband linear absorber 
(attenuator) of appropriate transmittance. As examples, 
we discussed CPA devices that perfectly absorb a chaotic 
laser signal and a frequency-modulated optical wave. 
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Appendix A: CPA for chaotic optical fields: 
derivation of the mean-field equations 

In this section we briefly derive the mean-field equa- 
tions (15) and (18) given in the text from the general 
maps (5) and (12) in case where the functional P de- 
scribes field propagation in a homogeneously-broadened 
two-level medium of length I. We will perform such 
a derivation for the more general case of the system 
S' of Fig. 1(b) with arbitrary injected signal /i(t) and 
transmittance Ta of the absorber; the laser system S of 
Fig. 1(a) is simply obtained after setting /i = 0, Ta = 0. 
Field propagation inside the ring cavity is described by 
the Maxwell-Bloch equations (see Sec. 4. 3 of Ref. [2]; see 
also Ref.dlj) 



dAi 

dz 
9Ai 

dni 



IdAi 



-gAi 



c dt 

-7_L [(1 + iA)Ai + mAi] 



(Al) 



-711 



n,-l-]^{AlK^ + A,K\) 



where: Ai{z,t), Ai(z,i), rii(z,t) are the intracavity elec- 
tric field envelope, polarization envelope and population 
inversion, respectively, normalized as in Refs.[2l [12]; 7|| 
and 7j_ are the population and dipole decay rates, re- 
spectively; g is the small-signal gain per unit length in 
the active medium; and A — (wo — ijj)/j± is the normal- 
ized detuning parameter between the cavity resonance 
frequency uj and the atomic transition frequency ujq. For 
the sake of definiteness, the gain mediunr (of lenght I) is 



placed at < z < I; similarly, we assume that the atten- 
uator Ta is a thin plate placed close to the beam splitter 
BS [see Fig. 1(b)]. In this way field attenuation after pas- 
sage across the absorber can be included in the scattering 
matrix of the beam splitter BS. This yields the following 
ring-cavity boundary condition for the field Ai{z,t) at 
the planes z — 0^ and z — C~ 

Ai (0+ , t) = Vrh (t) + VBTaAi {C-,t). ( A2) 

The output field Ei (t) escaping from the cavity through 
the beam splitter BS is then given by 



E,{t) = J^A,{Q\t)-^-^. 



(A3) 



Obviously ni and Ai in Eqs.(Al) are defined in the range 
< z < /, and for z > I (i.e. in the empty cavity region) 
the first equation of (Al) is still valid provided that we 
assume 5 = 0. The field Ai{£'~ , t) aX z — C~ and time t 
entering in Eq.(A2) can be computed as a functional of 
the field Ai (0+ , i - Tr) at plane z = 0+ and time t-Tj^, 
where Tr — C/c is the cavity transit time. In fact, after 
introduction of the new variables 
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= z , ri = t- z/c 


s.(Al) take the form 
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(A5) 



Equations (A5) can be integrated in the interval < 
£, < I and for < 77 < Tr once the initial conditions 
Ai{£_ = 0~^,r]) and Ai(^,77 = 0), 711(^,77 = 0) are as- 
signed (initial-boundary value (Goursat) problem [l7]). 
Typically, for a medium initially at equilibrium in the 
absence of the field, one can take Ai(^,7/ = 0) = 
and ni{S,,ri = 0) = 1. In this way, one can compute 
A{^ = C-,rj) = A{^ = ?, 77) for < ry < Tr. In the phys- 
ical space-time variables {z,t), this means that we know 
the value of Ai{C~ ,t) in the time interval Tr < t < 2Tr. 
With such a solution, from Eq.(A2) one can then calcu- 
late Ai{0'^,t) is the time interval Tj^ < t < 2Tr, which 
serves as an initial condition to integrate Eqs.(A5) in 
the interval < ^ < / and for Tr < rj < 2Tr, (sec- 
ond round-trip) with the appropriate initial values of 



^i{^:V = Tr) and ni{^,r] = Tr) computed at the pre- 
vious step. Iteration of such a procedure enables one to 
calculate ^1(0+,^) — ipiit) at any time t > 0, and thus 
the output field Ei{t) according to Eq.(A3). 

An important case is the one obtained by taking the 
limit gl — > 0. In this case, from Eqs.(A5) it follows that 
the variables ^i, Ai and rii are almost constant functions 
with respect to the ^ variable, i.e. ^1(^,77) ~ Ai(0, 7y), 
Ai(C:'7) "^ ^1(0,77) and ni{^,i]) ^ 711(0,77). The small 
change Ai{^ — C^,ri) — Ai{£^ — 0^,7;) can be computed 
from the first of Eqs.(A5) as Ai{£_ = £^,7;) - Ai(^ = 
0+,7;) ~ —glAi{^ = 0,77), which in terms of the space- 
time physical variables (z, t) reads 



AiiC~,t + Tr) ~ A,{0+,t) - glAi{0,t). 



(A6) 



Indicating by ?Ai(t) = Ai{z = 0+,t), Ai(i) = Ai(z = 0,t) 
and 77i(i) — ni{z — 0,t), from Eqs. (A2) and (A6) it 
follows that ipi (t) satisfies the delayed equation 

Mt+Tn) = VTf^it+TR)+^RTa[Mt)~9lMt)] (A7) 

which should be associated to the differential equations 
for Ai(t) and ni{t) 
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(A8) 



In their present form, Eqs.(A7) and (A8) represent a sys- 
tem of coupled differential-delayed equations that should 
be integrated with the initial conditions 77i(0), Ai(0) and 
ipiit) for < i < Tr. The common mean-field and single- 
mode laser model [H [2] is obtained from EciS.(A7) and 
(A8) by taking the further limits T = l-R^ 0,Ta ^ 1, 
with C = gl/T finite, and assunung that the injected and 
intracavity fields /i (t) and '0i (t) vary slowly with respect 
to time over one cavity transit time Tr. Under such as- 
sumptions, the delayed equation (A7) can be replaced by 
the following differential equation 



Mt) + TR 



dipi 
dt 



T 



VThit)+^aMt)-^Mt) 



-glKi{t) 
(A9) 



I.e. 



^ ~ -k[V'i(0 + 2CAi(i)] - Ki7/-i(t) + J^^ (AlO) 



where k and ki are defined by Eqs. (16) and (19) given in 
the text. The system of coupled differential equations 
(AS) and (AlO) is precisely the modified single-mode 
Lorenz-Haken laser model (18) given in the text. Note 
that the usual Lorenz-Haken laser model (15) is simply 
obtained after setting fi—Q and Tq = 1 (i.e. ki =0). 
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